Is subdiffusional transport slower than normal? 
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We consider anomalous non-Markovian transport of Brownian particles in viscoelastic fluid-like 
media with very large but finite macroscopic viscosity under the influence of a constant force field 
F. The viscoelastic properties of the medium are characterized by a power-law viscoelastic mem- 
ory kernel which ultra slow decays in time on the time scale r of strong viscoelastic correlations. 
The subdiffusive transport regime emerges transiently for i < r. However, the transport becomes 
asymptotically normal for i>r. ft is shown that even though transiently the mean displacement 
and the variance both scale sublinearly, i.e. anomalously slow, in time, (5x(t)} oc Ft a , (Sx 2 (t)) oc t a , 
< a < 1, the mean displacement at each instant of time is nevertheless always larger than one 
obtained for normal transport in a purely viscous medium with the same macroscopic viscosity ob- 
tained in the Markovian approximation. This can have profound implications for the subdiffusive 
transport in biological cells as the notion of "ultra-slowness" can be misleading in the context of 
anomalous diffusion-limited transport and reaction processes occurring on nano- and mesoscales. 

PACS numbers: 05.40.-a, 05.10.Gg, 87.16.Uv 



INTRODUCTION 



The widespread occurrence of anomalously slow diffu- 
sion and transport |l|-|9( in biological cells is still not com- 
monly appreciated in spite of a growing experimental evi- 
dence and support [lCj-|l8[ . One of the main psychological 
obstacles on the way to a wider recognition is that ultra- 
slowness seems intuitively be rather obstructive for the 
corresponding diffusion-limited biochemical reactions, es- 
pecially if it is caused by divergent mean residence times 
(MRTs) in trapping domains [l|, 0, 0, @] created by ran- 
dom meshwork of cell's cytoskeleton. From this perspec- 
tive, the occurrence of subdiffusion might be more asso- 
ciated with physics of dying, rather than with physics of 
life. Moreover, the bulk of experimental biophysical data 
is traditionally interpreted in terms of normal diffusion 
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and even modern biophysics textbooks (see e.g. 
discuss only normal diffusion, despite appreciating the 
fact of existence and importance of intracellular molecu- 
lar crowding which clearly obstructs diffusion via increase 
of the effective medium's viscosity 22|. The increase of 
effective viscosity depends also on the size of diffusing 
particles. So, "for molecules smaller than 1 nm, it's sim- 
ilar to that of water; for particles of diameter 6 nm ( such 
as a protein of mass 10 5 g/mol), it's about 3 times that of 
water. For 50-500 nm particles, it's 30-300 times that of 
water" 21] (p. 571). In this respect, typical globular pro- 
teins are in the range of 2-20 nm (diameter) [21, 22 1 and 
mRNA molecules are about 400-800 nm in diameter [h| . 
However, the traditional thinking and prejudges can also 
be the reasons for overlooking anomalous diffusion and 
transport regimes (probably mostly transient) as recent 
experimental work uncovers [13 Bill El- Moreover, 
the occurrence of subdiffusion clearly depends not only 
on the size of single macromolecules, but also on their 
concentration, i.e. on the degree of molecular crowding 
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One of approaches to anomalously slow diffusion and 
transport is traditionally based on the assumption of 
divergent MRTs in trapping domains [2, 0, UQ- Of 
course, MRT (r) in any finite spatial domain with lin- 
ear size Ax can never diverge in real life. However, it 
can largely exceed a characteristic diffusion time, rr> ~ 
(L 2 / Kq) 1 /" , required to subdiffusionally explore on av- 
erage, (Sx 2 (t)} oc K a t a with < a < 1, a finite volume 
with linear size L 3> Ax, where the correspond- 

ing subdiffusion coefficient measured in cm 2 /see" [23j |. 
Then, the approximation of infinite MRTs becomes phys- 
ically justified on the relevant mesoscopic scale L. For 
very large times t 3> (t), the diffusion becomes normal, 
(Sx 2 (t)) oc t. However, the corresponding spatial scale 
can largely exceed L and therefore the normal diffusion 
(and transport) regime can become of a little importance 
for certain processes in mesoscopic biochemical reactors 
of living cells, such as e. g. passive transport of mRNA 
macromolecules or large globular proteins [HI, fl(| , and 
in turn subdiffusion becomes of profound importance for 
such processes on mesoscopic scale. The approximation 
of divergent MRTs features the continuous time random 
walk (CTRW) approach to subdiffusion [H, 0, 0-0] . In this 
case, the position increments can be totally independent. 
Within the mean-field approximation, the CTRW trans- 
port is congruent [j] with jump-like transport in random 
potentials. Moreover, the overdamped continuous space 
Markovian Langevin dynamics in spatially varying po- 
tentials can be contracted onto such a semi-Markovian 
CTRW by doing properly spatial coarse-graining 3, 2|| . 
Then, a potential energy disorder can result in anoma- 
lous diffusion and transport in agreement with the semi- 
Markovian CTRW theory as recent work nicely demon- 
strates (see also paper by Lindenberg et al. 27 1 in 
this Special Issue). 

Alternatively, subdiffusion can result from the 
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medium's viscoelasticity [ 1 ll 1 1 3 . 16, 19, I28l - l3l| . It can be 
either due to viscoelasticity of the polymer actin mesh- 
work [n, m, m, m, or due to macromolecular crowd- 
ing in complex fluids 32[ as e.g. in cytoplasm of bacte- 
rial cells which are lacking static cytoskeleton 16,22, 3(| ■ 
Statistical analysis of the experimental single particle dif- 
fusional trajectories in bacteria [lj| reveals in fact the 
primarily viscoelastic origin of subdiffusion [33I ] . A main 
result of Ref. [HI is that the fractional Brownian motion 
scenario (see below) is more likely than one based on a 
semi-Markovian CTRW [34] . The authors of experimen- 
tal work [3l| came also to a similar conclusion. Within 
this alternative subdiffusional scenario all the moments 
of random time spent in finite spatial domains remain 
finite. The corresponding MRT is not only finite but it 
scales down to zero with Ax — > 0. The physical rea- 
son for subdiffusion here is very different. It occurs due 
to long-time anticorrelations in the position increments 
[30| . Considering Brownian particle of radius R, which 
starts to move at to = with velocity x(t) (we consider 
a one-dimensional case for simplicity), one expects it to 
experience a viscoelastic force 
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(t) = - / r](t - t')x(t')dt', 



(1) 



where 77(f) is a frictional memory kernel whose Laplace- 
transform fj(s) is related to the frequency-dependent 
medium's viscosity C(* w ) as vi s ) — 67ri?£(s). In the 
case of purely viscous fluids, and in neglecting the hy- 
drodynamic memory effects, £(t) = 2£q<5(£), where Co is 
the fluid's macroscopic viscosity, so that 77(f) = 2rj S(t), 
where 770 = 67ri?Co is the Stokes viscous friction coeffi- 
cient. For weakly viscoelastic fluids, £(t) = (o^exp(— vt), 
exponentially decays in time with rate v, and correspond- 
ingly 77(f) = 67ri?Co i/ exp(— vt) = «cxp(— v£), where k has 
dimension of a linear elastic force constant. This corre- 
sponds to the Maxwell theory of viscoelasticity (35| who 
derived the phenomenon of viscosity from medium's elas- 
ticity by assuming that the linear elastic force, F c \(t) = 
—n[x(t) — x(0)], acting on the particle can relax in time 
with rate v, yielding a viscoelastic force, i.e. F v _ e i(t) = 
—Kx(t) — vF v - e i(t). Indeed, if the force relaxation is 
very fast with respect to the change of particle's ve- 
locity, then F v _ e i(t) = — J Q r](t — t')x{t')dt' w —T]ox(t), 
with rjo = k/v, whereas in the opposite limit, F v _ e i(i) « 
—n[x(t) — x(Q)}. For strongly viscoelastic media one ex- 
pects that 77(f) decays in time much slower than expo- 
nential and a power law decay 77(f) oc t~ a can serve 
as a better model. For a fluid-like environment the ef- 
fective macroscopic viscosity Co = J °° ({t)dt should re- 
main, however, finite. It can be very large, but yet finite. 
Therefore, a long-time cutoff to the power-law must exist. 
In 1936, A. Gemant proposed a class of power-law vis- 
coelastic models which are consistent with this demand 
36l 37 1 . Its particular representative corresponds in the 
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f)(s) = / 77(f) exp(— st)dt = 
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in our notations. Here, 77(C)) — r\ is an effective asymp- 
totic friction coefficient and r presents a long-time mem- 
ory cutoff. The corresponding memory kernel is approx 
imately 
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for f <C t, where ?7 Q = r\r a ~ x , and T(x) is the familiar 
gamma-function. For i > t, 77(f) decays also in accor- 
dance with a power law, 77(f) oc f"~ 2 , i.e. elastic corre- 
lations are still rather strong. However, the correspond- 
ing integral converges ensuring that the asymptotic fric- 
tion coefficient 77 is finite. In the limit of infinitely large 
medium's viscosity yielding 77 — > 00, and infinitely long 
memory range, r — > 00, with ?7 a = 7/t q_1 kept constant, 
Eq. ([3]) becomes exact, fj(s) = 77 Q s Q_1 , and 
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where the last equality defines fractional Caputo deriva- 
tive of the order < a < 1 38] acting on x(t). r\ a 
can be named the fractional friction coefficient. Clearly, 
Eqs. (|3I4[) can serve as a good approximation only for the 
times <Ct. In the focus of this Letter is but the entire 
time evolution, interpolating between transient subdiffu- 
sion and asymptotically normal diffusion behavior. For 
example, r can correspond to the time scale of seconds or 
minutes, and then subdiffusion emerges on the time scale 
from microseconds to seconds or minutes, as in biological 

ceiis [IM3,[ii. 



SIMPLE MODEL 

We continue with a non-Markovian generalized 
Langevin equation (GLE) description 39r41| for an over- 
damped Brownian particle neglecting inertial effects. 
Then, 



v (t-t')x(t')dt' = f( Xl t)+at), 



(5) 



where f(x, f ) is a generally nonlinear force acting on the 
particle and £ (f ) is a thermal random force of the environ- 
ment. It is Gaussian, unbiased on average, and obeying 
the fluctuation-dissipation relation, 
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at the environmental temperature T. This is required for 
the consistency with thermodynamics at thermal equi- 
librium. In the above-mentioned limit r — ¥ 00, i] — > 
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FIG. 1. Mean displacement (in units of Fr/n) versus time (in 
units of r) under the influence of constant force F for several 
different values of a and the same rj, r. The limit of normal 
diffusion is achieved asymptotically from above. Transiently 
subdiffusing particles always cover larger distances than nor- 
mally diffusing particles with the same asymptotic frictional 
force constant 77. 



00, rj a — const, GLE ([5J is named also the fractional 
Langevin equation [9j upon the use of the corresponding 
abbreviation (jU for its lhs. In this limit, t;(t) is noth- 
ing else but the fractional Gaussian noise by Mandelbrot 
and van Ness 42| which presents an instance of 1// noise 
with the spectral power density S(u>) oc l/w 1_Q . Notice 
that generally the lower integration limit in Eq. ([5]) is 
t{) — > -co. It can be replaced, however, with to = since 
we assume that the particle starts to move at this time 
being initially localized, i.e. v(t) — x(t) = for t < t . 

Let us consider the transport under the influence of 
constant force F. Then, the above GLE can be easily 
solved formally using the Laplace-transform method, for 
any memory kernel rj(t). Transforming back to the time 
domain for the memory kernel ©, one obtains for the 
averaged mean displacement the simple result, 
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This exact solution is compared in FigJT]with the solution 
of the ordinary Langevin equation with memoryless fric- 
tion for the particle which experiences the same frictional 
force for the whole time span as our particle asymptoti- 
cally, or, said differently, the result of the Markovian ap- 
proximation to the considered dynamics. Clearly, for all 
times our particle moves in fact faster, covering larger 
distances and approaching gradually the limit of nor- 
mal diffusion from above. The instant time-dependent 
ensemble-averaged velocity (v(t)) := (Sx(t)) = Fjr\ + 
F/[n a T(a)t 1 ~ a ] is also always larger than its asymptotic 
value F/rj (a spurious singularity at t = can be elimi- 
nated, if to take the initial inertial effects into account). 
Mathematically, this is simply because x a 3> x, for x <C 1 
and < a < 1. The physics is also clear. In the limit 



77 oc Co — ► 00 the normally diffusing particle is get local- 
ized, (Sx(t)} — > 0, whereas our particle still moves, but 
ultra slow (per definition) since (Sx(t)} oc (t/r)". There- 
fore, the "ultra-slow" moving particle can cover larger 
distances. Furthermore, for any memory kernel in the 
studied model and for arbitrary constant F the variance 



of the particle position. 
391 



(Sx 2 (t)) = (x 2 (t))-(x(t)) 2 obeys 



<fa 2 (i)> = ^(Sx(t)), 



(8) 



and therefore it follows to the same pattern as in Eq. ([7} 
and Fig. [TJ 



(6x 2 (t)) = 2mt + 2K a t a /T{l + a) 



(9) 



Here, 

fact IS fractional diffusion coefficient related to tem- 
perature and fractional friction coefficient by the gen- 
eralized Einstein-Stokes relation, n a = ksT/rja, which 
contains the standard one, k,\ = ksT/n, as a particular 
case. 

We suppose that our observation is rather general. For 
example, the results in Ref. [H[ seem to agree with our 
line of reasoning. Indeed, mRNA macromolecules have 
in the related experiments radii in the range of 200-500 
nm. Furthermore, the normal diffusion coefficient in wa- 
ter was found to be K\ — I /jm 2 /sec (see Supplementary 
Material in [Hj]). From this, given the water viscosity 
Q w = 0.9- 10~ 3 Pa-sec, one can estimate the corresponding 
radius as R = ksT/ (6tt( w ki) which gives R « 242 nm for 
T = 300 K. Let us assume that R sa 250 nm. Then, the 
corresponding macroscopic normal diffusion coefficient in 
cytosol should be by the factor of r ~ 300 smaller than 
one in water [2lJ (see the above quotation in Introduc- 
tion). This yields K{ yt) w K X /r « 3.33 • 10~ 3 /jm 2 /sec. 
However, the experiment yields not normal but subdif- 
fusion with a « 0.7 (see Fig. 2(a) in [1511 and n a in 
the range from 10~ 3 to 10 -2 /im 2 /sec 0,7 [34{. Assuming 
K a = 10 -2 /j,m 2 /sec 7 for this value of R (smaller parti- 
cles in experiment should also subdiffuse faster) one can 
conclude that subdiffusion can indeed cover larger dis- 
tances than normal diffusion with jr. Further- 
more, one can estimate the transition time r. Given the 
relation r = (r//?7 Q ) 1/(1 - Q) = (kJk^) 1 ^ 1 -^ which 
follows within our model one obtains for it t « 55 sec. 
This is a rather reasonable estimate since subdiffusion 
regime lasts in those experiments up to 30 sec, cf. Fig. 
2(a) in [Hj]. 



CONCLUSIONS 

The discussed phenomenon might seem paradoxical, 
even though its explanation is almost trivial. Neverthe- 
less, it has profound implications for subdiffusion in bi- 
ological cells. First of all, the occurrence of subdiffusion 
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on some transient time scale r and the corresponding 
mesoscopic spatial scale L ~ (2n a T a ) 1 ^ 2 does not con- 
tradict to the bulk of macroscopic experimental data in- 
dicating typically a normal diffusion 21|. Even more 



important, the overall transport is in fact faster than its 
long time normal asymptotics that results from a drastic 
reduction of the effective normal diffusion coefficient for 
large macromolecules due to molecular crowding effects 
in cytosol. Very important is also the fact that on the 
time scale t <C r, by use of the memory kernel (j4} as 
approximation, the resulting random process is nothing 
else the fractional Brownian motion (fBm) [3^, 42]. The 



fractal dimensionality of the fBm trajectories is known 
to be df — 2/ a [HI, i.e. df = 2 for the normal Brown- 
ian motion. However, for a < 2/3 and for 3d Euclidean 
embedding dimension it becomes df =3 (a fractal tra- 
jectory cannot have larger Hausdorff dimension that the 
Euclidean space in which it is living). This means that 
a subdiffusionally searching particle explores much more 
thoroughly finite volumes than a normally diffusing par- 
ticle on the same spatial and time scales. Therefore, 
transient subdiffusion can provide only advantages for 
the diffusion-limited reactions in tiny biochemical reac- 
tors of living cells, which are densely stuffed with differ- 
ent macromolecules, without some principal drawbacks. 
This is rather unexpected and paradoxical conclusion. 
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